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The feasibility of detecting a heavy charged Higgs boson, mjj± > rrit + rn-^, decaying in the 
H^ — > tb channel is studied with the fast simulation of the atlas detector. We study the 



(^ ' 99 ~^ H^tft production process at the lhc which together with the aforementioned decay chan- 

^— ^ \ nel leads to four 6-quarks in the final state. The whole production and decay chain reads 

f^ I gg — > f6H^ -^ tibb -^ bbbblvqq' . Combinatorial background is a major difficulty in this multi-jet 

f^ ■ environment but can be overcome by employing multivariate techniques in the event reconstruc- 



tion. Requiring four 6-tagged jets in the event helps to effectively suppress the Standard Model 



Q^' backgrounds but leads to no significant improvement in the discovery potential compared to anal- 

i^ ' yses requiring only three 6-tagged jets. This study indicates that charged Higgs bosons can be 



discovered at the lhc up to high masses (?7iH± > 400 GeV) in the case of large tan/3. 
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I. INTRODUCTION 

The only particle predicted by the Standard Model (SM) that has so far not been de- 
tected is the Higgs boson. In the Minimal Supersymmetric extension to the Standard Model 
(MSSM) [l|, 12] the Higgs sector is enlarged to contain 5 particles: 3 neutral (h°,H°)A°) and 
two charged (H+,H^) Higgs bosons. Whereas the detection of one neutral Higgs bosons 
would be compatible with both the SM and the MSSM, the detection of a charged spin-0 
particle such as the charged Higgs boson predicted by the MSSM would unequivocally point 
towards new physics beyond the SM. This note describes the potential of the atlas exper- 
iment to detect a heavy charged Higgs boson, i.e. a charged Higgs boson heavier than the 
top quark, decaying in the H^ -^ th channel. 

Other experiments have also searched for the charged Higgs boson and set lower limits 
on the charged Higgs boson mass. The combined LEP experiments provide a preliminary 
exclusion of charged Higgs bosons with mH± < 78.6 GeV at the 95 % CL '3| . At the Tevatron, 
CDF and DO searched for the charged Higgs boson in the decay of top quarks produced 
in pp — > ti reactions. These searches exclude the low and high tan/3 regions up to charged 
Higgs masses of ^ 160 GeV y, Isf- 

The Higgs sector of the MSSM is determined by two free parameters at tree level, most 
often chosen to be the mass of the CP-odd neutral Higgs boson, rriA, and the ratio of the 
vacuum expectation values of the two electroweak Higgs doublets, parametrised by tan/3. 
The decay modes of the charged Higgs boson in the MSSM are given as a function of the 
charged Higgs mass in Figure ^ for two different values of tan/?, 1.5 and 30. These plots 
show that the main decay channel of heavy charged Higgs bosons for r?7,H± > mt + mb is the 
decay into a top- and a 6-quark . However, searches in this decay channel have to resolve 
the problem of a large multi-jet background. For this reason the most promising channel for 
the search for the charged Higgs boson heavier than the top quark is the H+ -^ tVj. decay 
channel, as it provides a lower background environment p. 

The H^ -^ th decay channel, assuming mH± ^ ''^t + '^b, has been studied in a previous 
note y, using the 2 — > 2 production process gh — * H^'^t and detecting 3 6-jets in the final 
state. However, as was shown in that report, the large background from Standard Model 
tt-production complicates the detection of the charged Higgs boson and limits the mass 
region for a charged Higgs discovery to masses below about 400 GeV for low or high values 
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FIG. 1: Branching ratios for different decay channels of the charged Higgs hoson in the MSSM 
as a function of the charged Higgs mass. The upper plot assumes tan/3 = 1.5, and the lower plot 
tan/? = 30. 

of tan/3. The purpose of the present study is to try to extend the discovery reach beyond 
this hmit. 



Recently it was suggested [8| that by utihzing the 2^3 production process gg -^ H^tb 
in combination with the H"^ -^ th decay, the fourth 6-jet inherent in the signal process could 
be detected, resulting in a greater rejection of the Standard Model background processes. 
We therefore study a heavy charged Higgs boson in the production and decay chain 
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thR^ -^ tthh -^ bbbbluqq', 



(1) 



where one of the top quarks is required to decay leptonically in order to provide a hard 
isolated lepton to trigger on. The SM background processes that lead to the same final state 
with four 6-tagged jets are 

99 -^ tibb (2) 
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FIG. 2: (T X BR of the process gg — > tbU.^ —>■ ttbb —>■ bbbbluqq' for a charged Higgs mass of 
mH± = 300 GeV as a function of tan/3. 



and 



99 -^ ttgg + ttqq, 



(3) 



where, in the latter case, two of the hght jets are misidentified as 6-jets. 

Analyses searching for a charged Higgs boson in the production processes gb -^ H^t and 
gg -^ H'^tb generally suffer from a lack of sensitivity for intermediate values of tan/5. The 
relevant part of the MSSM Lagrangian describing the H^t6 Yukawa coupling is given by 



C 



-\/2mw± sin 9w 



H+ (mbtan/3 tbR + rrit cot/3 1 bi) 



(4) 



which has a minimum at tan/? = Jmt/m\y. This behavior is illustrated in Figure El showing 
the cross section times branching ratio (BR) for process (^ as a function of tan/?. The dip 
around tan/5 = Jmi/m\^ ^ 7 is apparent. 

In the following section the event generation and detector simulation are described. Sec- 
tion mil describes the analysis which is divided into two likelihood selections presented in 
the sections fill Bl and llll (^[ respectively. Section HVl summarises the results, and conclusions 
and an outlook are given in section ^\ 



II. EVENT GENERATION AND SIMULATION 

The event generation and detector simulation for the various Monte Carlo (MC) samples 
used in this analysis are done within the athena framework in the atlas Software Re- 
lease 6.0.3. 

The signal process (^ is generated for a charged Higgs boson mass range of 
mH± = 200 - 800 GeV using HERWIG 6.5 jsl llo|. Table HI hsts the mass points at which 
MC samples are produced and gives production cross sections for a set of tan/5 values. In 
the cross section calculations the factorisation and renormalisation scales, fip and fin, re- 
spectively, are set to the mean transverse mass so that fip = fij^ = J2i=b,t{PT + "^D/2 = 
{m^{t) +m'^{b))/2 = {m"^). As all cross sections are calculated at leading order and no next 
to leading order calculations exist to date, an optimal choice of the QCD scale is not obvious. 
The choice of this particular scale is guided by the demand to use comparable scales in the 



signal and background calculations for consistency and by results obtained in [11[ although 
the process considered there can not be compared directly to the signal process considered 
in this report. Choosing the QCD scale is one of the main systematic uncertainties when 
predicting the signal and background cross sections. The value assumed here provides esti- 
mates of the cross sections to be expected but can by no means be considered as definitive. 
Other choices of the QCD scale or mh evaluation may result in cross sections differing by up 
to a factor 2. This topic is further discussed in section HVl 

The strong coupling constant ag is evaluated at the 1-loop level and a running b- 
quark mass is used. A central value of rrit = 175 GeV is assumed for the top quark 
mass and the CTEq5l parton density function is used throughout the analysis. H+ —>■ tb 
branching ratios are evaluated with HDECAY 3.0 1^ where the decay to supersym- 
metric particles is switched off (0FF-SUSY=1). We evaluate the branching ratios 
in the Maximal Mixing Scenario as described in |13j, assuming the top quark mass 
mentioned above: MgusY = 1 TeV, Ms = 200 GeV, fi = -200 GeV, Mg = 800 GeV and 
At = Ab = Xt + fi/tanP with Xt = -\/6Msusy- The branching ratio of the W^ decaying 
to quarks is assumed to be 2/3, and the BR to a lepton (electron or muon) plus the ac- 
companying neutrino is taken to be 2/9. Since each of the W^ can decay leptonically or 
hadronically, an overall factor of 2 has to be applied. This leads to the following relation: 
a X BR = a^ncl. x BR(H+ -^ tb) x 8/27. 



TABLE I: Summary of signal MC samples for the various charged Higgs masses studied in this 
analysis. Inclusive and exclusive cross sections are quoted for a selection of tan/J values, assuming 
fj!p = fij^ = (?n,^(t) + m'j.{b))/2 and a running b~quark mass. Each MC sample contains 10^ gen- 
erated events. 



mH± 


tan/3 


CTind. [Pb] 


cr X BR [pb] 


200 
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0.17 
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1.3 


0.31 


300 
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0.35 
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0.35 
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20 


0.11 


0.029 


500 


30 


0.13 


0.034 


600 


40 


0.13 


0.033 


800 


30 


0.024 


0.0063 



TABLE II: Summary of background MC samples studied in this analysis. 



Process 




Generator 


CTind. [pb] 


generated events 


gg -^ ttbb 




AcerMC 1.2 


10.3 


10 M 


qq -^ ttbb 




AcerMC 1.2 


0.61 


IM 


gg ^ Z/j/W - 


-> ttbb 


AcerMC 1.2 


1.1 


IM 


tt + jets 




HERWIG 6.5 


405.0 


50 M 



Samples for the background processes gg — > ttbb, qq — > ttbb and gg -^ Z/'-f/W —>■ ttbb are 
produced with AcerMC 1.2 [l4'| in stand-alone mode. The QCD energy scale is chosen 
such that Qqcd = ^(Pt + "^i)/4 = {rn^). These samples are then passed to HERWiG 6.5 
within the Athena framework for fragmentation and hadronisation. In order to study the 
systematic uncertainty due to different fragmentation schemes, the gg -^ ttbb sample is also 
passed to pythia 6.203 jl^ for a similar treatment. The effect of mis-tagging light jets 
as 6-jets is studied with the help of a large tt +jets sample generated with herwig 6.5. 
Table |n] summarises the background samples and their inclusive cross sections. 

The ATLAS detector is simulated with the fast detector simulation atlfast as it is rep- 



resented in the atlas Software Release 6.0.3. This package is based on the FORTRAN 
implementation of the same package [l6|. Jets are reconstructed with a cone based algo- 
rithm using a cone size of AR = 0.4. Only jets having a minimum transverse momentum 
of pt > 10 GeV and lying in the pseudorapidity range of \ri\ < 5 are accepted for this anal- 
ysis. An efficiency of 90 % to identify isolated charged leptons is assumed. Jet energy and 
momentum calibration and the 6-tagging of jets is performed within the atlfastb routine 
of the ATLFAST simulation package. The possibility to tag a jet as a 6-jet is limited by the 
inner tracker acceptance range of |?7| < 2.5. A 6-tagging efficiency of 60 % is assumed when 
simulating samples for the low luminosity option of the LHC, 50 % for the high luminosity op- 
tion. Rejection factors of Re = 10 and Rj = 100 are chosen for c- and light jets respectively. 
The 6-tagging efficiencies and rejection factors are static, i.e. they do not depend on the 
pseudorapidity t] or transverse momentum px of the jets. All plots and tables shown in this 
analysis refer to the low luminosity option of the LHC, assuming an integrated luminosity of 
30 fb~^ unless explicitly stated otherwise. 

A. Jet Part on Matching 

In order to construct and test the performance of the event reconstruction algorithm (see 



section IIIIB|1 . it is necessary to know the link between a generated parton and a detected 
jet or lepton. The former is often referred to as the "Monte Carlo truth", and the latter will 
be referred to as "reconstructed objects" in the following. Initially no such link between a 
parton and a reconstructed object is provided by the MC generator or the detector simulation 
program and the association is far from straightforward. In this analysis the problem is 
handled approximately by solving the assignment problem as described in [13] ■ The quantity 
which is minimised is the sum of all distances between the generated partons after final 
state radiation (FSR) and their associated reconstructed object 4-vectors. The distance 
between the 4-vector of a parton and the 4-vector of a reconstructed object is given by 
AR = y/A'^T] + A^0, the distance in pseudorapidity-azimuthal angle space. If the distance 
between a parton and its reconstructed object exceeds 0.4 it is assumed that no association 
is possible. Further, no association is attempted if any of the initial quarks after FSR has a 
4-momentum outside the acceptance range of the JetMaker algorithm in atlfast. 



III. ANALYSIS 

The analysis has three parts. In the first step all events are required to pass a set of 
cuts in order to reject most of the SM background and to ensure the minimum prerequisites 
needed for subsequent reconstruction. 

The second part is intended to find the combination of jets that correctly reconstructs 
the two top quarks and the charged Higgs present in the final state of the signal process. 
For each event the most likely correct combination is found with the help of a selection 
procedure described in section IIIIBI This likelihood is referred to as the "combinatorial 
likelihood". 

Once the correct combination is found for each event, a second likelihood selection, the 
"selection likelihood" described in section IlllCt aims at separating the signal from the SM 
background processes. 

A. Preselection 

In the preselection, events with a topology clearly distinct from the signal topology are 
rejected. This ensures that only the main backgrounds discussed in section H] need to be 
studied further. The preselection requires: 

• exactly 1 isolated lepton (/ = e or /i) with transverse momentum p^ > 25 GeV, 
Pj. > 20 GeV and pseudorapidity \ri\ < 2.5, 

• exactly 4 5-jets with \t]\ < 2.5 and pt > 20 GeV and 

• at least 2 light jets with \ri\ < 5 and px > 20 GeV. 

In order to trigger on the signal events the detection of a high-py lepton is required. The 
cuts applied to the p^ and t] of the isolated lepton are chosen such that they meet the 
requirements of the ATLAS trigger system. When running in the high luminosity option the 
cut on the jets' transverse momenta is increased from pj^ > 20 GeV to pj. > 30 GeV. The 
efficiency of the precuts for signal events depends on the assumed charged Higgs mass and 
ranges from 1.78 % for mH± = 200 GeV to 4.41 % for mH± = 800 GeV. The precut efficiencies 
are summarised in Table ITlll 



TABLE III: Preselection efficiencies for the different charged Higgs masses studied in this analysis. 
The efficiencies are quoted for the low luminosity (LL) and the high luminosity (HL) option of the 

LHC. 

mH± [GeV] preselection efficiency, LL [%] preselection efficiency, HL [%] 

200 L783 ± 0.013 0.5687 ± 0.0075 

250 3.171 ±0.018 1.121 ±0.011 

300 3.577 ±0.019 1.281 ±0.011 

350 3.738 ±0.019 1.436 ±0.012 

400 3.932 ±0.019 1.509 ±0.012 

500 4.123 ±0.020 1.644 ±0.013 

600 4.281 ±0.020 1.715 ±0.013 

700 4.363 ± 0.020 (not studied) 

800 4.411 ±0.021 1.872 ±0.014 

In order to reconstruct the leptonically decaying W^ (Wj^p) the 4-monientum of the 
daughter neutrino needs to be reconstructed. The x- and y-components of the neutrino 
momentum are assumed to coincide with the measured missing transverse momentum com- 
ponents p^^^^ and p™'^^ respectively. The z-component however can not be measured but 
must be calculated by solving the equation: 



m 



^± = (E^ + Elf - (p; ± p^)^ with E^ = \pl\. 



This equation can result in two or zero real solutions for p^. If two solutions are found, both 
are kept for later evaluation in the event reconstruction algorithm. However, in approxi- 
mately 25 % of the events no solution is found. In order to keep these events and still be 
able to reconstruct the leptonically decaying Ww in those otherwise fatal cases, the collinear 
approximation approach described in 1^ is adopted: p^ = pf is assumed if no solution can 
be found, and the resulting W^ 4-momentum is rescaled to match mw± • This increases the 
Ww reconstruction efficiency from 75 % to 100 % and only a small loss in the resolution of 
the reconstructed leptonically decaying tiep is observed. 



B. The Combinatorial Likelihood 

The final state of tlie signal process (HJ is quite complex, featuring four 6-jets, two light 
jets from the hadronically decaying W^(Wj|3^j) and an isolated lepton plus missing transverse 
momentum from the leptonically decaying W^(Ww). Quark and gluon jets from initial and 
final state radiation and the underlying event are also present, increasing the jet multiplicity. 
Initially it is unknown which reconstructed objects should be combined to reconstruct the 
two W^s, the two top quarks, and finally the charged Higgs boson. The combinatorial 
likelihood aims at identifying the correct reconstructed objects to combine and thereby 
making the correct reconstruction of the whole event possible. In order to incorporate as 
much information available from each event as possible a multivariate technique is chosen to 
find the correct combination of reconstructed objects for each event. We choose to implement 
a likelihood selection distinguishing two classes where the first class represents the correct 
combination and the second class all the wrong ones. The likelihood formalism used in this 
analysis is outlined briefly in the following, generalising to n classes of events: 

For each of the m observables Xj used to distinguish between the n classes, the normalised 
probability density functions (pdf ) 

fi{xi), where i = 1, . . . ,m and j = 1, . . . ,n (5) 

have to be determined. The probability that an event belongs to class j when measuring 
the value Xi for variable i is given by 

The likelihood C that an event belongs to class j when measuring m variables xi, . . . , Xm is 
then given by the normalised product of the probabilities Eq.lHlover all m variables: 

m 

a = -^^ (7) 

k=ll=l 

The likelihood has values in the range 0,. . . ,1. Information about possible correlations 
between the input variables is neglected by this procedure. 

When constructing the pdfs for the combinatorial likelihood it is essential to know the 
correct association of partons to reconstructed objects as described in section IIIAI The 
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FIG. 3: One of the possible Feynman diagrams for the signal process gg -^ tbH^ -^ bhbblvqq' , 
illustrating the labelling of the partons adopted in this analysis. 

fraction of signal events passing the precuts for which a vahd association to partons is found 
depends on the charged Higgs mass and ranges from 45 % at 77iH± = 200 GeV to 58 % at 
mH± = 800 GeV. Only events having a valid parton association can be used to obtain the 
pdfs for the correct- and the wrong-combination class. 

Any algorithm used to reconstruct the events has a chance to find the correct combina- 
tion only if the correct four jets are 6-tagged and the two light jets originating from the 
hadronically decaying W^ pass the precut constraints. Of those events passing the precuts 
and obtaining a valid association to partons, only approximately 65 % fulfill this require- 
ment. This means that for ^ 35% oi signal events passing the precuts the completely correct 
reconstruction is doomed from the beginning. 

The combinatorial likelihood is based on the following 9 variables, where the labelling of 
partons in the signal process is illustrated in Figure El 

1. rrijj: the invariant mass of any two light jets. For the correct combination this mass 
should be within the W^ mass peak around ■m-\!i^± = 80.4 GeV whereas the distribution 
of invariant masses of pairs of jets not originating from a W^ is rather flat. 

2. rrijji,: the invariant mass of any two light jets and one of the four 6-jets. The correct 
combination should reproduce the top mass. This variable aims at correctly recon- 
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structing the hadronically decaying top quark thad- 

3. rriwb'- the invariant mass of the isolated lepton, one solution for the neutrino recon- 
structed as described in llll Al and one of the 6-jets. This variable aims at reconstructing 
the leptonically decaying top quark tiep. 

4. pT{p2): the pr of the &-jet assumed to originate from the charged Higgs decay. 

5. pribo): the pr of the assumed companion 6-jet produced in the gg -^ t6H^ process. 

6. AR(j, j): AR between any two light jets. Like rrijj this variable helps to reconstruct 



W^ 



had" 



7. AR{jj,b): AR between the sum of any two light jets and a 6-jet. Like rrijjb this 
variable is related to the reconstruction of t^ad- 

8. AR(Z, b): AR between the isolated lepton and a 6-jet. Like mi,jh this variable aims at 
reconstructing tiep. 

9. AR.{b2,ti): AR between the 6-jet and the top quark candidate originating from the 
charged Higgs decay. For the top quark candidate all possible modes of reconstructing 
a top quark are considered. 

The probability density functions for the nine variables used in the combinatorial likelihood 
are shown in Figure ID for a charged Higgs mass of 400 GeV. Overflow bins are included in 
the normalisation of the histograms and used when calculating the combinatorial likelihood 
output. The corresponding normalised likelihood distributions for the correct-combination 
class are shown in Figure Elfor the correct combination and all the wrong ones. As expected, 
the distribution corresponding to the correct combination peaks at 1 whereas the distribution 
representing all the wrong combinations peaks at 0. However, it is important to note the 
tail in the distribution representing the wrong combinations up to high likelihood values. 
The total number of combinations of the reconstructed objects to reconstruct the event 

completely is given by 

(m \ 

where m is the number of light jets in the event and N^, is the number of solutions for 
the neutrino. The 4! represents the number of possibilities to order the four 6-jets and 
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FIG. 4: The prohability density functions for the nine variables used in the combinatorial likelihood 
for a charged Higgs mass of mjji = 400 GeV. All distributions are normalised to unity including 
overEow bins. Correct combinations are represented by a solid line, all wrong combinations by a 
dashed line. 

the factor of 2 reflects the possible associations of Wj^p and W{|!j^ to the top quarks. This 
number N depends on the number of light jets in the event and is generally quite large. 
Hence the number of wrong combinations is large and the possibility of one of those wrong 
combinations having a combinatorial likelihood output higher than the correct combination 
is not negligible, thus reducing the probability of identifying the correct combination. 
For each event the combination yielding the highest correct-combination likelihood is 
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FIG. 5: The normalised combinatorial likelihood output distributions for the correct-combination 
class. The correct (solid) and the wrong (dashed) combinations are shown for a charged Higgs 
mass of rrajji = 400 GeV. Only combinations with a likelihood output larger than 0.7 are accepted. 

treated as the correct combination. Nevertheless, if this selected combination yields a likeli- 
hood value below 0.7 the event is rejected. The efficiency of this cut varies between 90 % for 
mH± = 200 GeV and 95 % for mH± = 800 GeV for the signal process and is approximately 
85 % for the main gg — * tibb background. 

The performance of the combinatorial likelihood is checked using the Monte Carlo truth 
information to associate the final state partons with the reconstructed objects as described 
in section Hi Al An event is classified as correctly reconstructed if the four 6-jets and the two 
light jets are correctly associated with their corresponding final state partons and the correct 
lepton is found to be isolated. Some performance benchmarks of the combinatorial likelihood 
are shown in Figure IHl as a function of the charged Higgs boson mass. The squares indicate 
the fraction of correctly reconstructed hadronically decaying W^, referred to as purity in the 
following. This purity does not depend strongly on the charged Higgs mass and lies between 
53 % and 60 %. The purity of reconstructing the two top quarks is represented by the open 
circles. Here the purity depends strongly on ?7iH± and rises to values above 25 % only for 
charged Higgs masses either very close to rrit or for ?nH± > 700 GeV. Similar behaviour can 



14 



^ 



u 



100 
90 
80 
70 
60 
50 
40 
30 
20 
10 



I 










■ hadronic W 




I 










O top-quarks 
A whole event 
)K charged Higgs 


z 










I 














"■■" > ^ ^ ^ 


z * 










= 
A 






*. 
* 


E 






\ ^ ^ c> 

* i *■ 


E 


« * « * A ^ 
A A A A 




~ 1 1 1 1 




1 1 1 1 1 1 1 1 1 



100 



200 



300 



400 



500 



600 



700 



800 



900 



nijj. [GeV] 



FIG. 6: Purity of the combinatorial likelihood in reconstructing the hadronically decaying W^ 
(squares), the two top quarks (circles), the charged Higgs boson (stars) and the whole event (tri- 
angles) over the considered charged Higgs mass range. 

be seen for the purity of the charged Higgs boson (stars) and the whole event reconstruction 
(triangles). The only variables used in the combinatorial likelihood that depend strongly 
on mH± are variables 4 and 9. The p^ of the 6-jet from the charged Higgs decay depends 
strongly on the mass of the charged Higgs boson and is well separated from the average Pt 
of the other 6-jets in the event only for very small or quite large charged Higgs masses. A 
similar statement applies to the distance AK(b2,ti) between the top quark and the 6-quark 
originating from the charged Higgs decay. A light charged Higgs boson is produced with 
a sizable boost and its decay products will have a small distance in AR-space. On the 
other hand, a heavy charged Higgs boson is produced nearly at rest and hence the distance 
between its decay products will be large. The mass dependence of the performance of the 
combinatorial likelihood is mainly determined by the mass dependence of variables 4 and 9. 
It should be noted that in the mass region mH± = 250 - 600 GeV the purity of the charged 
Higgs reconstruction does not exceed 30 %. As a consequence the reconstructed charged 
Higgs mass is substantially blurred by the combinatorial background. Hence the detection 
of a clear mass peak in the reconstructed charged Higgs mass distribution is difficult. The 
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FIG. 7: Reconstructed masses for the hadronically decaying W^ and the hadronic and leptonic top 
quark for events passing the cut on the combinatorial hkelihood. The signal is generated assuming 
mH± =400 GeV. 

correct reconstruction of the whole event is important for the performance of the selection 
likelihood discussed below, since it relies on the correct association of reconstructed objects 
to the charged Higgs boson and the top quarks. The ability to distinguish signal from back- 
ground processes is already diminished by the imperfect performance of the combinatorial 
likelihood. 

Before combining a W^ and a 6-quark to form a reconstructed top quark, the W^ 4- 
momentum is scaled to reproduce mw± = 80.4 GeV. The reconstructed hadronic W^ and 
the two reconstructed top masses for events passing the cut of 0.7 on the combinatorial 
likelihood output are shown in Figure [7| Reconstructed charged Higgs masses are shown for 
mH± = 200, 400 and 800 GeV in Figure |H1 The solid line represents the reconstructed mass 
obtained with the combinatorial likelihood. To predict the detector performance and to 
illustrate the effect of the combinatorial background the same distributions are also shown 
using the Monte Carlo truth information to select the correct combination of reconstructed 
objects as the dashed lines. The effect of the combinatorial background is clearly visible as 
a tail in the reconstructed charged Higgs mass distribution, especially toward higher recon- 
structed masses. For higher charged Higgs masses detector effects become more prominent. 
Even when the MC truth information is included a large tail toward lower reconstructed 
masses develops. 
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FIG. 8: Reconstructed charged Higgs masses for ?tt-h± = 200, 400 and 800 GeV. The sohd hne shows 
the reconstructed mass as obtained with the combinatorial lil^elihood. The dashed line illustrates 
the charged Higgs mass resolution obtained when utilizing the Monte Carlo truth information to 
match the tree-level partons to the reconstructed objects. 

C. The Selection Likelihood 

To enhance the signal and suppress the Standard Model background a second likelihood 
selection is implemented. The selection likelihood distinguishes three classes of events: 1) the 
gg -^ tbE."^ signal process, 2) the gg/qq — > tibb background, and 3) the gg -^ Z/j/W -^ tibb 
background process and is implemented using the same formalism as described in sec- 
tion IIIIBI It exploits differences between the distributions of the signal and the Standard 
Model backgrounds in the following four variables: 

1- '"^feob2- ^he invariant mass of the two 6-jets not originating from a top quark decay. In 
the signal events one of the jets originates from a heavy charged Higgs boson whereas 
in the background processes both 6-jets originate from gluon splitting. Hence the 
invariant mass of the two jets is expected to be lower in background than in signal 
events. 

2. cos^(6o, ^2): the cosine of the angle between the two 6-jets not originating from a top 
quark decay. Since the two 6-jets originate from gluon splitting for the background 
processes they are expected to be collinear whereas the distribution should be flat for 
the signal process. 

3. cos6'(6o + ^2): the cosine of the azimuthal angle of the 60 + &2 jet system. 

4. cos6'(tboost, H^recon): the cosine of the angle between the reconstructed charged Higgs 
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boson momentum and the reconstructed top quark associated with its decay, where 
the reconstructed top quark 4-momentum is boosted into the charged Higgs boson 
rest frame. 

The corresponding probabihty density distributions are shown in FigurelHlfor ??ih±=600 GeV. 

Variables related to transverse momenta or invariant masses of jet systems tend to shift 
the background peak in the reconstructed charged Higgs mass distribution towards the 
nominal charged Higgs mass. Therefore only variables involving angular correlations are used 
in the selection likelihood. The only exception is m^gbj for which it has been demonstrated 
that no such shift occurs for cuts of less than 0.4 on the resulting selection likelihood. 

The normalised selection likelihood output distributions for the signal class and assuming 
mH± = 600 GeV are shown in Figure^! These distributions are calculated for each charged 
Higgs mass under consideration. Signal events are separated from the Standard Model back- 
ground by selecting only those events yielding a selection likelihood output larger than 0.2. 
This optimal cut value is found by varying the cut on the likelihood output in steps of 0.05 
and requiring the method to yield the highest discovery potential over the selected range of 
charged Higgs masses. Figure ITT] shows the resulting reconstructed charged Higgs mass dis- 
tribution for a choice of charged Higgs masses. Here an integrated luminosity of £ = 30 fb~^ 
and tan/3 = 80 is assumed. Only a slight shift of the peak in the reconstructed charged Higgs 
boson mass for the backgrounds is observed for growing mH± • For '^h* ~ 400 GeV the peaks 
in the reconstructed charged Higgs mass for the signal and the background processes can be 
separated. 

All events within a mass window of ±100 GeV around the nominal charged Higgs mass 
are selected. The width of the mass window has little influence on the discovery potential 
and hence is not optimised. The number of selected signal and background events is then 
treated like in a simple counting experiment and the Poisson significance is calculated for 
each charged Higgs mass and value of tan/5. The resulting 5 a discovery contour is presented 
and discussed in the next section. 

Figure ^1 summarises the signal selection efficiency and the expected number of back- 
ground events after all cuts have been made for the whole range of charged Higgs masses 
studied. The signal selection efficiency lies within 1.4 — 2.7% and reaches its maximum 
around mH± ~ 300 GeV. 
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FIG. 9: The probability density functions for the four variables used in the selection likelihood 
for a charged Higgs mass of r?T,jj± = 600 GeV. All distributions are normalised to unity. For each 
variable the distributions corresponding to the gg -^ tbH signal process (solid), the gg/qq -^ tibb 
background process (dashed) and the electroweak gg -^ Z/j/W -^ tibb process (dotted) are shown. 

IV. RESULTS 

The results of the analysis are described in this section in terms of 5 o" discovery contours 
in the {rriA, tan/3) plane. They are presented for integrated luminosities of 30 fb~^ for the low 
luminosity option and 300 fb~ for the high luminosity option of the LHC. In the latter case a 
6-tag efficiency of e^ = 0.5 is assumed and the p™° cut on all jets is raised to pj- > 30 GeV. 
The degradation of jet-energy measurements due to pile-up is taken approximately into ac- 
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FIG. 10: The normalised selection likelihood output distributions for the signal class. Distri- 
butions for the gg —>■ tbH.^ signal process (solid), the gg/qq — > tibb background process (dashed) 
and the electroweak gg -^ Z/^/W — > ttbb process (dotted) are shown for a charged Higgs mass of 
mjji = 600 GeV. Only events yielding a selection likelihood output larger than 0.2 are selected. 

count by choosing the high luminosity option in the atlfast simulation package. Figure IT^ 
shows the expected discovery contours taking no systematic uncertainties into account. The 
charged Higgs boson can be detected for tan/5 values down to 35 for mH± ~ 250 GeV based 
on an integrated luminosity of 30 fb~^. For the high luminosity option and 300 fb~^ the reach 
in tan/5 goes down to approximately 28 for the same charged Higgs mass region. The anal- 
ysis presented here depends heavily on the 6-tagging performance and the reconstruction of 
jets with a relatively low px. This explains why only a small improvement in the discovery 
potential is observed when switching to the high luminosity option. 

The uncertainty in the prediction of the signal and background cross sections due to the 
choice of QCD scale and running or pole 6-quark masses is quite large as already discussed 
in section IHI Figure ^^ shows some cross section predictions obtained with herwig 6.5 
for the gg -^ tbH^ process, assuming different QCD scales and nib evaluations. The lower 
three curves represent cross section predictions for a running 6-quark mass and renormal- 
isation and factorisation scales of /i^ = /i|; = {jriTit) + mr(6))^/4, (m'^(t) + m|-(6))/2, and 
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FIG. 11: Reconstructed charged Higgs masses for rrajji = 200, 400, 600 and 800 GeV and tan/3 = 
80. The histograms are normahsed to the expected event rate for C = 30fb~ . The gg -^ t6H^ 
signal process (yellow) is shown on top of the Standard Model backgrounds gg/qq — > ttbb (green), 
gg -^ Z/-f/W -^ ttbb (light blue) and tt + jets (pink). 

(mH± + rrit + rnbY' . Wliereas only a small difference is observed between the predictions of 
the first two choices, a rather large reduction in the expected signal cross section is observed 
if a QCD scale of rmji + mb + mt is assumed. However, NLO calculations for the 2 — > 2 

^^ ri 



process gh -^ tH 



19| and the gg -^ ttH process |11] show that this choice of scale might 



be too high. The same studies show also that the cross sections are likely to be overesti- 
mated when using a pole 6-quark mass. We therefore adopt a running 6-mass, ensuring 
K-factors larger than 1. To illustrate the effect of a larger signal cross section prediction 
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FIG. 12: The signal selection efEciency and the total expected number of background events after 
all cuts in the mass window of ±100 GeV assuming an integrated luminosity of 30 fb^ . The dots 
represent the charged Higgs masses that were studied, and the line is a smooth interpolation. 

we nevertheless show the cross sections expected when assuming a 6-quark pole mass in the 
uppermost curve in Figure El and the corresponding improvement in the discovery contour 
in Figure IT^ f left plot). The latter plot shows that improvements in the discovery potential 
due to K-factors > 1 might be sizable. 

The main qqlqq —>■ tibb background cross section prediction is also very sensitive to the 
QCD scale [ij] and the uncertainties on the cross section prediction are of the same order as 
for the signal process. However, here we will assume that the background cross section can 
be measured using side-bands in the reconstructed mass distribution which are relatively 
signal-free. The precision of this procedure depends on the charged Higgs mass and on 
the integrated luminosity available. No detailed study is conducted here, but to give some 
indication of how the discovery potential is affected by this uncertainty on the expected 
Standard Model background, we assume an uncertainty of 5-10 % in the background nor- 
malisation, guided by the studies done in [l8|. If the background MC samples produced 
with AcerMC are passed to pythia for string fragmentation and hadronisation instead 
of herwig's cluster fragmentation, differences between 5 % and 10 % are observed in the 
background prediction, depending on the charged Higgs mass. 

To illustrate the effects of uncertainties on the Standard Model background prediction of 
this order of magnitude, we show the effect of 5 % and 10 % uncertainties on the discovery 
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FIG. 13: The 5 a discovery contours for the process gg —>^ tb}i^ ^ tibb ^ bbbbluqq' in the 
(mA) tan/3) plane for the low and the high luminosity option of the lhc, assuming integrated 
luminosities of 30 fb~ and 300 fb~ respectively. A common renormalisation and factorisation 
scale of fi'jp = fj,^ = (?7i^(t) + m^(&))/2 and a running b-quark mass are assumed when evaluating 
the cross sections of the various processes involved. No systematic uncertainties are taken into 
account. 

potential in Figure E| Again, the corrections are found to be sizable. 

V. CONCLUSION AND OUTLOOK 

This note analyses the discovery potential for a charged Higgs boson heavier than the top 
quark produced in the 2 — > 3 process gg -^ t6H^. The subsequent decay to heavy quarks 
H+ -^ tb is considered, leading to a final state consisting of four 6-jets, two light jets and 
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FIG. 14: cross section predictions for the gg -^ t6H process obtained with herwig 6.5 assuming 
different QCD scales and b-quark mass evaluations. A value of tan/3 = 30 is chosen. 
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FIG. 15: Influences of systematic uncertainties on the discovery potential assuming an integrated 
luminosity of 30 fb^ and fx'p = fij^ = {mj,{t) + mj.{h))/2. The left plot shows the improvement in 
discovery potential if a pole b-quark mass is assumed instead of a running mi, when evaluating the 
gg -^ t6H^ cross section. The right plot shows 5 a discovery contours taking systematic uncertain- 
ties on the background cross section normalisation of 5% and 10 % into account. Here a running 
mb is used. 
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one electron or niuon plus missing energy. The whole production and decay chain reads 
gg -^ thYi"^ -^ tibb -^ bbbbluqq'. Studying the 2^3 process offers the possibility to detect 
four 6-jets in the final state and thereby reduce the Standard Model background considerably 
compared to case in which the 2^2 production process gb -^ tH^ is considered. 

One of the main difficulties to overcome when trying to reconstruct signal events is the 
high number of possible combinations of paired reconstructed objects in order to reconstruct 
the charged Higgs boson. It is shown that the reconstruction is possible by employing multi- 
variate techniques, in which angular correlations are also taken into account. By employing a 
likelihood selection separating the Standard Model background from signal events, it is pos- 
sible to detect the charged Higgs boson in the gg —^ tbH^ — > tibb -^ bbbblvqq' channel down 
to values of tan/5 ^ 28 for charged Higgs masses around 250 GeV assuming an integrated 
luminosity of 300 fb~^. However, the theoretical uncertainties related to the cross section 
predictions for both the signal process and the main Standard Model background are quite 
large and lead to a sizable uncertainty in the expected discovery contour in the (r/iA, tan/3) 
plane. NLO corrections to the signal cross section might result in an improvement in the 
discovery potential whereas expected uncertainties when measuring the Standard Model 
background contribution will degrade the result. We present the 5 a discovery contour in 
Figure IT^ using a running 6-quark mass and ^\ = ^\ = {m'^{t) + m\{b))/2 and taking no 
systematic uncertainties on the background normalisation into account. 

The goal of this analysis was to utilise the detection of the fourth 6-jet in the signal 
process in order to extend the discovery region for the charged Higgs boson at high charged 
Higgs masses as suggested in js| . This analysis shows that the encouraging results obtained 
in [8| do not hold when detector effects and mis-tagging of 6-jets are more properly taken 
into account. 

A direct comparison to a previous analysis J7| where the 2 — i> 2 production process 
gb — > tYii^ was used to produce a heavy charged Higgs boson which subsequently also decays 
to heavy quarks, H+ -^ tb, is not possible at this stage, since the cross section predictions 
and production mechanisms for the Standard Model backgrounds that are assumed in the 
two cases are different. 

Finally it should be noted that the results presented here might be subject to another large 
systematic effect. As was mentioned in section |Hj the 6-tag efficiencies and rejection factors 
assumed are static, i.e. they do not depend on r] nor on the pt of the jet under consideration. 
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This is clearly a rather crude approximation especially in the present analysis for which the 
detection of 4 6-jets is crucial. More reliable results should be possible in the future using 
a more accurate (pr,'7)^dependent parametrisation for the 6-tagging efficiency. 
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